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Abstract: While the computed transported turbulence dissipation rate, s, works well as part of a differential 
equation-based turbulence model in predicting turbulent flows, it doesn’t seem to work well when used to 
determine the Kolmogorov length-scale (l K oi) which, like the other Kolmogorov scales, exists within the viscous 
sublayer portion of the inner turbulent boundary layer zone. Using smay lead to an increase in l Koi as roughness 
increases, the opposite of what should happen. It is proposed here to replace the computed s (and its level at the 
1 st point off the wall as dictated by wall function 5 ) Tvith\ the bne resulting from basic law-of-the-wall sublayer 
relationships which includes the Prandtl-Schlidtiting (P-S‘) roughness effect. This approach enables physically 
correct prediction ofb Koi , particularly a reliable decrease th&rebfwith increasing roughness level. 
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I. Introduction 


The Kolmogorov length-scale (l Ko i), under high Re conditions, is commonly expressed as (v 3 /£) ly/4 but 
the computed decay rate of turbulence kinetic energy, 8, doesif’t seem to be adequately sensitive to surface 
roughness, leading to dVer-estimation angUometimes even wrong-behavior Of l Koi !This could be a consequence 
of both the 8 transport equation and its wall boundary condition lacking correct roughness effects, even though si 
is set by the law-of-the-wMl including roughness effect. 

An alternative way^toj^Xpress r 4 KjDi is proposed here, its basic formulation being desensitized to whether the first 
off-wall computational points f&yfwithin the viscous sublayer (including the solve-to-wall level of y + < 1) or in the 
logarithmic overlap B (see Appendix for flow property details in these regions).The basic aim is to select formulation 
ingredients such that the resulting l Koi is minimal (to admit more small roughness levels) and sensitive to the 
roughness size since.suriace (boundary) conditions are strongly influencin|, the immediately adjacent viscous 
sublayer where Kolmogorov micro-scales "are defined. Effect of roughness height on boundary layer features is 


seen below: 


\ 



Figure 1. Velocity profile vs. roughness height 


As seen in Fig. 1, the larger roughness element, r 2 , causes increase in t w (red broken line velocity profile, U 2 ) but 
the outer velocity portion of U 2 is the same as that of the smaller roughness, ri. Conservation on mass-flow-rate, 
m, dictates adjustments within the lower section of the boundary layer. Thus the inner portion of the boundary 
layer (viscous sublayer and logarithmic overlap) shrinks in the normal-to-wall direction, whence the dissipative 
eddies become smaller. And since l Koi is a measure of these eddies, it is expected to also become smaller as the 
roughness increases. 

Following the next Theory section, five examples are given to exhibit the usefulness of the proposed approach to 
compute l KoX . 
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Theory 

1. Definition of l Koi 
Following Eq. (Cl) in the Appendix, 
2 A £ k 2 

8 vR y {l - e~ AeR y) 
where 


A £ = C^ /4 /(2k) 

R y =min{y + /C 4/4 ,^y +2 }, 4* s 0.0266 

Here the 2 nd term within the parentheses represents the viscous sublayer. Thus 


l Kol 


£ J 


, 1/4 

\ -Z. 

R y (l - e A £ R yJ 

) Vfc 

2A £ 


1/4 


( 1 ) 


(la) 

(lb) 


( 2 ) 


where k (see Appendix Eq. (Al) and Fig. 2) is chosen at its maximum value, attained within the logarithmic layer: 


fc 



(3) 


Importantly, x w magnitude increases with roughness size, thereby increasing k and, by Eq. 2, decreasing l Koi . This 

/ v 3\ 1/4 

is concluded also for the original l Koi = (—) based on the behavior of 8 within the viscous sublayer, leading 

to t Koi = (y) = (2i4 fe ) -1 / 4 v J[~") ■ Fig- 3 shows experimental data (Bellucci et al., 2014) of increasing k 

with roughness level in the viscous sublayer. Since £~ k 2 /v within the sublayer, s also gets larger with increasing 
roughness. Consequently, both original and proposed t KoX will be reduced in size as roughness increases. This has 

(Ky) l/ 4 y3/4 

been found also in Marati et al., 2004, where l Koi = —- jm. 

(jw/P w ) 

The effect of roughness on t w is based on the P-S (Prandtl and Schlichting, 1934) skin friction formula when the 
1 st computational points are within the viscous sublayer: 

C f = [2.87 +l.S8log 10 ^) /(l+ v -^Ml) ,r w - -\ Qx UlC f and u T = V(r/p) w (4) 


This is combined with equations 2 and 3, thus introducing roughness effect explicitly. 



2. The Kolmogorov length-scale as a function of s 
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Based on the Appendix Eqs. (B7-B8), the Kolmogorov length-scale grows along the surface streamline distance, 
s, as follows: 


Kolsublayer 


(pA)w — 


yjv w {du/dy\ 



V 

1 

N. 

* 

■layer 

_V ( J /<?)w_ 

L(du/dy) w J 


1/4 


(5a) 

(5b) 


and since (du/dy) w commonly becomes smaller along s, the Kolmogorov length-scale grows along it as seen in 
Fig. 4. 


y 



(du/dy) w 

Figure 4. (9U t /5y) w diminishes along the streamwise direction 


Test cases 

The commercial flow solver CFD++ (see details in Chakravarthy, 1999) was used with the realizable k-s closure 
(Goldberg et al., 1998) to compute the following flow examples, using a non-equilibrium wall function (see 
Appendix) which includes the effects of dp/ds and roughness, the latter by Eq. (4). All grids used for the following 
flow cases have been previously established to yield mesh-independent solutions. 


(I) Flat plate 

Air flow over a rough plate was computed with two roughness levels: Ks=3 micro (p)-meter and 0.6 
mm. l st off-plate grid centroids were once within the viscous sublayer (y*=7) and once in the 
logarithmic layer (y*=56). Fig. 5 compares t Koi based on the standard, directly computed 8, 
(v 3 /^) 1 / 4 ,with that based on Equations 1-4. The choice of Ks = 3xl0' 6 m emanates from Schlichting 
1968, p.612, where camouflage paint of this equivalent sand-grain height was applied to a war plane, 
resulting in a substantial increase in friction drag relative to that of the unpainted airplane. The 
standard Kolmogorov length-scale is about 5 times larger than this Ks, rendering it negligible even 
though in reality it has a large effect on viscous drag. This indicates that the standard t KoX formula 
may not be adequate for cases involving rough surfaces, suggesting the alternative formulation 
hereby proposed. The standard formula predicts a considerably larger drop in t Koi than the proposed 
method does as the Ks magnitude increases to 0.6 mm (solid vs. dashed lines in Fig. 5F). A similar 
trend is observed in the y + =56 case (Fig. 5R), however, since the Kolmogorov scales exist only 
within the viscous sublayer, the 1 st computational point must be inside this sublayer, therefore y + =56 
is not acceptable. Since the smaller /^is predicted by the current proposal, it admits a larger range 
of roughness sizes. This test case indicates that limiting the 1 st off-surface computational points 

within the viscous sublayer (including y* < 1 for direct solve-to-wall) is the preferable approach 
(Fig. 5F) which adheres to the physics since the Kolmogorov scales exist only within the viscous 
sublayer. 
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Kolmogorov length-scale Kolmogorov length-scale 

Flat plate, M = 0,125, y :+: = 7 Flat plate, M = 0,125, y + = y :+: = 56 




Figure 5. l Koi profiles along a rough plate, Ks=3pm and 0.6 mm. 1 st computational points being inside the viscous 

sublayer(L) and within the logarithmic layer (R) 

(ii) Flow over a curved surface 

This is a case of flow over a curved surface (from Suga et al., 2006). Wall geometry and computational mesh are 
shown in Fig. 6(L), consisting of 4672 cells of which 2426 are triangles and the rest are quadrilaterals adjacent to 
the wall. 91 cells are along the wall and 50 normal to it, of which 23 are the quadrilaterals. The 1 st row along the 
wall has a 0.5 mm height, corresponding to y*= 10 upstream of the curve and y*= 8 along it. The wall is subject 
to air flow at a temperature of 288 K and Uoo=22 m/s. 

Figure 6(R) shows l Koi profiles along the wall for rough surfaces with Ks=0.003 mm and 0.6 mm. As in the flat 
plate case, Eq. (2) predicts lower l Koi levels than those from the standard approach, thereby widening the range 
of admissible small Ks levels. Both methods predict lower levels of l Koi under the higher roughness level as 
expected. 

Kolmogorov length-scale 


convex wall, U = 22 m/s, y* = 8 - 10 




Figure 6. (L) Computational domain, (R) /^profiles 


It is noted in both the above examples (as well as in the next ones) that the leading-edge zone of the surface 
exhibits very large l KoX as predicted by the standard method, t Koi = (v 3 /^) 1 / 4 , likely due to underprediction of 
c in that region, whereas the proposed approach (Equations 1-4) behaves considerably better in this zone. 


(Ill) NACA 0012 airfoil at a=6° 


Manuscript id. 752206930 


www.iistre.com 


Page 4 
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Addressed in Lu and Liou, 2009, this airfoil mesh has 113362 triangles and 12680 quadrilaterals. There are 947 
quadrilateral cells around the airfoil with y*=5 in the trailing edge region on both pressure and suction sides and 
y*= 15 in the leading-edge zone (Fig. 7L). 

Kolmogorov length-scale 

NACA 0012 U = 232 ft/s, 6°, y* = 5 




Figure 7. (L) mesh around the airfoil, (R) l KoX profiles along the NACA 0012 suction surface. 


As in previous cases, the proposed formulation yields a smaller /^than the standard one does and in this case 
also produces no oscillation at the leading-edge zone (Fig. 7R), as opposed to the standard approach. Comparing 
the 0.6 mm rough wall results with those of the 3 p-meter one shows the consistency of the proposed approach 
versus the inconsistency of the standard one regarding the behavior of £ KoL when switching from lower to higher 
roughness: Along the trailing edge section l Kol is predicted larger for the higher roughness whereas the opposite 
is physically correct. On the other hand, the proposed method (red lines) is consistently correct along the entire 
airfoil. 


(iv) Water flow over a sand dune 

This test case is a sand dune with water flowing over it (from Suga et al., 2006). The flow includes a separation 
zone at the bottom of the dune. Equivalent sand-grain roughness heights of 3 p-meter and 0.6 mm are compared 
as before. The domain is 1.6 m long and the free air/water surface is located 0.292 m above the level of the 
beginning of the dune, treated as a symmetry boundary to prevent flow normal to it. The bulk mean velocity is U 
= 0.633 m/s and the Reynolds number, based on U and the free surface height, is 175,000. The quadrilaterals grid 
size is 16,121 and the 1 st computational points above the dune are within the viscous sublayer, y*=7 being a typical 
level. Figure 8(L) shows the grid and 8(R) shows plots of /^profiles along the dune for the two roughness levels, 
as predicted by the standard and proposed formulations. 

Kolmogorov length-scale 

water over sand dune, U = 0.633 m/s, y :+: = 7 




Figure 8. (L) Sand dune mesh, (R) l Koi profiles along dune (lower surface) 

It is noticed in Fig. 8R that the proposed l Koi shows no difference in spite of a factor of 200 between the two 
roughness levels. The reason is the high density of water. For example, at the mid-point (x=0.8 m) the two levels 
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of roughness give rise (Eq. 3) to A k = = 0 031 = 3 x 10 5 (m/s) 2 and that imposes a negligible change in 

Pvp 1003.9 ' 

t Koi , Eq. 2. However, like in previous examples, the new l KoX is smaller than the one from the standard approach 
and it doesn’t suffer from over-predicting l Koi in the inlet zone as the standard method does. 

(V) Channel flow with a lower rough wall including a curved portion 

Here is a channel flow with a rough lower surface which includes a curved portion. Song and Eaton, 2002, describe 
this flow case. Figures 9 show topology and grid for this case. Fig. 9(L) includes details like the circular curved 
ramp section length of L=70 mm, its radius of 127 mm and its 21 mm height. The mesh is especially dense in the 
curved ramp area due to the flow being separated there. The upper channel wall is smooth while the lower wall is 
rough (originally with Ks=1.2 mm but here Ks=3 p-meter and 0.6 mm are used). Inflow speed is 20 m/s and 
turbulence levels are Tu=0.5% and p t /p=30. The 11,200 cells mesh has 8,000 quadrilaterals and the rest are 
triangles. The 0.15 mm height of the next-to-walls grid layer translates into 7> y*> 4 on the ramp wall. 

Fig. 10 compares l Koi along the lower surface for Ks=0.6 mm and Ks=3 p-meter levels of roughness, using both 
standard and new methods. Both approaches reduce t Koi as the roughness increases, with the standard method 
producing a considerably larger jump as in previous cases. However, the proposed method predicts a smaller l KoX 
level for the 3 p-meter roughness than the standard one predicts for the 0.6 mm roughness (a trend seen in previous 
flow examples) and, as before, it avoids the standard method’s trend of overpredicting i Ko {\n the inlet region. 



TO 



Figure 9. (F) Topology and dimensions, (R) Grid 


Kolmogorov length-scale 
emved rairup, U =. 20 in/s-, y* - 4 .7 



Figure 10. /^profiles along curved ramp (lower surface) 


Concluding observations 

(1) The larger the roughness, the thinner the inner boundary layer - including the sublayer, due to increased 
(<9U t /dy) w (Fig. 1). That translates into smaller dissipative eddies. Since l KoX is a measure of these eddies within 
the sublayer, it is also expected to shrink in size with increasing roughness. The proposed approach adheres to 
this: l KoX decreases when moving from a low roughness level to a higher one (examples 1, 2, 3 and 5). However, 
the standard method predicts the opposite trend in the trailing edge region of example 3 (airfoil), this being 
unphysical. On the other hand, in examples 1, 2, 4 and 5 the standard method does predict the correct trend. 
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(2) When switching from Ksi= 3 |Li-meter to Ks 2 = 0.6 mm roughness, a change in f Ko iof — 

Zkol, 2 

/2.87+1.58iO£ 10 ^V' 25 

(-) is expected, based on Eqs (2-4). With s=l m this ratio is about 1.6. In thel st example the 

y2.87+1.58 log 10 -^— J 

standard method yields 1.8 and the new one gives 1.2. In the 2 nd one the standard method predicts 1.6 while the 
proposed one yields again 1.2. In the 5 th case 2.0 results from the standard approach whereas the new one predicts 
a ratio «1. Thus, the standard method is closer to the expected ratio, however, it exhibits a wrong-oriented ratio 
in example 3. 

(3) Again, in the 3 rd example (airfoil), along the trailing edge zone, the standard method predicts that the rougher 
wall causes a large increase in l Koi whereas the opposite is correct, based on fundamental near-wall relationships 
from the law-of-the-wall in the viscous sublayer (Equations 1-4 and the Appendix). Along the rest of the airfoil 
the standard method does recover the correct behavior, but the proposed approach maintains correct prediction 
trends along the entire airfoil. 

(4) The standard approach introduces high levels of l KoX in the inflow (or leading edge) region, as seen in cases 
1, 2, 4 and 5. The new method doesn’t suffer from this phenomenon due to avoiding usage of the solved variable 


(5) The proposed method is based on the fundamental rules of the sublayer’s wall function (Equations 1-3 which 
already include roughness effects by using the P-S [Prandtl and Schlichting, 1934] skin friction formula, Eq. 4) 
and those rules have been proven correct since many years ago, based on numerous experimental data. Since there 
is no such thing as a truly smooth surface, wall functions have long been valid for mildly rough walls as long as 
roughness is taken into account in the form of a sand-grain equivalent height, without changing the surface 
geometry. Of course, large roughness (e.g. Ks > 3 mm) should be implemented as part of the surface geometry 
and grid, in which case the currently proposed treatment isn’t necessary. 

(6) Whenever there are doubts, it is always advantageous to resort to fundamental relationships. In the present 
case eliminating 8, computed from its own transport equation and subject to its own wall BC (with £ 1 = 2A k u\ /v, 
see Appendix B, and Pk=0 in the viscous sublayer) and replacing 8 with the fundamental wall function-based 
description thereof (Equations 1-3 which include the P-S roughness-influenced Cf, Eq. 4) is a safe and secure way 
of dealing with l Kol . 


Nomenclature 


A k =0.0266, 
B = 5.5, 

C f 


Q=0.09, 

E=9.8, 

Ks 

k 


M 

Pk 

P 

Re 

Rt 

Ry 

s 

s 

Tu 
U, u 
U t 

u + =U t /u T , 

u T =(kl/p) w 1/2 , 

y 

y + =yU T /v, 


constant used in the viscous sublayer 
constant in the logarithmic wall function branch 
skin friction 

coefficient of eddy viscosity in some k-s turbulence models 

constant in some forms of the logarithmic Law-of-the-Wall 

near-wall damping function 

sand grain height representing roughness level 

turbulence kinetic energy 

length-scale 

Mach number 

production of turbulence kinetic energy 
static pressure 
flow Reynolds number 
turbulence Reynolds number 
turbulence Reynolds number based on y + 
strain rate 

distance along body contour 
freestream turbulence intensity 
streamwise velocity component 
tangent-to-wall velocity component 
non-dimensional speed for wall functions 
friction speed 
normal-to-surface distance 

normal-to-wall dimensionless coordinate for wall functions 


y*=C^ A s[ky/v, k-based normal-to-wall dimensionless coordinate for wall functions 
Greekjetters 

s dissipation rate of turbulence production 
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y =Cp/C v , 

specific heat ratio 

k^0.41, 

von Karman constant 

F 

dynamic molecular viscosity 

Ft 

dynamic eddy viscosity 

< 

ii 

p 

kinematic molecular viscosity 

Vt=p t /p, 

kinematic eddy viscosity 

p 

density 

T 

subscripts 

shear stress 

1 

denoting 1 st off-wall computational point 

00 

freestream conditions 

c 

corrected for compressibility effects 

Kol 

of Kolmogorov 

t 

turbulent or tangential 

w 

at wall 
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Appendix: BASIC SUBLAYER/LOG-LAYER RELATIONSHIPS 


A. Logarithmic layer 


k = 

w t /'J V (r/p)w 


(Al) 

£ = 

u|/(Ky) 


(A2) 

V t = 

Kyu r 


(A3) 

5 = 

u T /(Ky) 


(A4) 

Kolmogorov velocity scale: (/]g) 

II 

(A5) 

/c 

Realizable time-scale: T t = - - 

£ 

* ? 

II 

(A6) 

Ry ~ 

= //c 


(A7) 

R t = 

■■ C~ V \ Ry 


(A8) 
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u + = ln(e B y +1/K ^, B = 5.5, k = 0.41 

(A9) 

The above assumes that f^=l in the expression v t = C^f^k 2 / e. 

B. Viscous sublayer 
u + = y + 

(Bl) 

s-m"- 4,2 00266 

(B2) 

k = A k (u T y + ) 2 = A k u 2 

(B3) 

2 A k u* 

£ = - 

V 

(B4) 

v t = C^A^vy + \ A^ = 0.0084 

V2 

(B5) 

Ur 

s = — 

V 

(B6) 

Kolmogorov velocity scale: (ye) 4 = (2A k )^sfvS = 0.48VvS = 

1 

^(7^=0.341/^ (B7) 

1 v 

Realizable time-scale: T t = ? 

y/AkU r 

(B8) 

CM 

+ 

II 

aT 

(B9) 

D _ di v + 4 — 3 

H t ~ 2 y ~ 2 

(B10) 

C. General e formula (based on the above) 


k 2 

£ ~ 2A£ vR y a-e- A * R y ) 

(Cl) 

cT 

(C2) 

R y = mm{C” 1/4 y + , V^fcy +Z } 

(C3) 

or 


R y = max{ A /2R t; 2i4 £ R t ) 

(C4) 
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